UNCLASSIFIED 


AD  NUMBER 


AD822698 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution:  Further  dissemination  only 
as  directed  by  Navy  Weather  Research 
Facility,  Naval  Air  Station,  VA  23511,  Jun 
1965  or  higher  DoD  authority. 


AUTHORITY 


USNPS  ltr,  19  Sep  1972 


THIS  PAGE  IS  UNCLASSIFIED 


AD822698 


NWRF 


SYNOPTIC  ANALYSIS  AND 
FORECASTING  OF 
SURFACE  CURRENTS 

A  REPRINT  OF 

FLEET  NUMERICAL  WEATHER  FACILITY 
TECHNICAL  NOTE  NO.  9 


NAVY  WEATHER  RESEARCH  FACILITY 

BUILDING  R-48,  NAVAL  AIR  STATION 
NORFOLK,  VIRGINIA-23511 


SYNOPTIC  ANALYSIS  AND  FORECASTING 
OF  SURFACE  CURRENTS 


Fleet  Numerical  Weather  Facility 
Monterey,  California 


Technical  Note  No.  9 
June  1985 


By 

CAPT  W.  E.  Hubert,  USN 
and 

T.  Laevastu 


FOREWORD 


"Synoptic  Analysis  and  Forecasting  of  Surface  Currents"  as  reproduced 
herein  is  an  edited  version  of  the  paper  distributed  previously  in  limited 
quantity  as  Fleet  Numerical  Weather  Facility  Technical  Note  No.  9. 

This  publication  presents  operational  techniques  for  the  analysis  and 
forecasting  of  surface  water  currents  from  synoptic  data.  The  difficult  fore¬ 
cast  problem  has  been  reduced  to  manageable  proportions  by  treating  the  total 
current  as  several  distinct  components  at  different  time  and  space  scales,  and 
by  describing  the  individual  prediction  relationships  ir.  terms  of  basic  con¬ 
siderations  and  procedures. 
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INTRODUCTION 


In  navlg&tlonj  fisheries  and  other  sciences  connected  with  the  sea.  It  la 
often  necessary  to  forecast  or  hlndcast  the  direction  and  speed  of  surface  cur¬ 
rents  at  a  specified  time  and  place.  Past  attempts  at  current  prediction  have 
had  only  moderate  success.  It  can  he  argued  that  progress  on  the  quantitative 
understanding  sr.d  prediction  of  currents  has  beer,  hampered  by  a  tendency  to 
consider  ocean  currents  to  be  analogous  to  vind  currents  In  the  atmosphere. 
Actually,  It  appears  that  ocean  currents  are  affected  by  an  even  greater  num¬ 
ber  of  factors,  and  that  different  techniques  may  be  required  for  successful 
prediction.  This  summary  attempts  to  avoid  past  shortcomings  by  treating  the 
total  current  In  terms  of  a  number  of  components.  The  existing  quantitative 
knowledge  of  cause-effect  relations  of  these  components  is  first  reviewed,  and 
procedures  are  suggested  for  computation  of  currents  in  different  time  and 
space  scale->. 

A  summary  of  techniques  for  surface  current  prediction  has  been  distilled 
from  a  voluminous  literature  search.  Qnphaals  will  be  put  on  the  basic  consid¬ 
erations  and  procedures  rather  than  attempting  a  mathematical  reformulation  of 
the  problem.  The  theoretical  treatises  on  large  scale  (oceanic  to  global)  cur¬ 
rent  depiction  will  not  be  discussed  here;  good  reference?  to  thes*  works  can 
be  found  in  textbooks  and  other  summaries  (e.g.,  Robinson  [95]).  Descriptions 
of  some  of  the  major  ocean  currents  can  be  found  in  books  by  Stocncl  [109],  and 
Reid  [94]  . 


This  sunnary  has  been  built  upon  an  earlier  draft  by  one  of  the  authors; 
however.  It  has  been  extensively  revised  on  the  basis  of  recent  experiences 
in  analysis  and  forecasting  of  surface  currents  and  their  verification  at  the 
U.  S.  Fleet  Numerical  Weather  Facility. 


1.  GENERAL  NATURE  AND  VARIABILITY  OF  SURFACE  CURRENTS; 

PRESENT  STAGE  OF  KNOWLEDGE  AND  PREDICTABILITY 

The  literature  on  ocean  currents  contains  a  number  of  seemingly  contra¬ 
dictory  descriptions  and  theories.  Many  of  the  controversies  arise  from  the 
differences  which  are  observed  In  current  behavior  (both  in  space  and  time). 
These  observed  differences  result  in  most  cases  from  the  various  combinations 
of  driving  or  influencing  forces  which  prevail.  For  this  reason,  it  is  be¬ 
lieved  that  ocean  currents  can  best  be  understood  and  characterized  by  de¬ 
scribing  their  space/timo  variability  through  study  of  their  principal  com¬ 
ponents  . 

Currents  are  typically  unsteady  in  direction  and  speed  (even  in  cases 
where  the  tidal  components  are  subtracted  and  where  the  local  winds  vary  lit¬ 
tle).  T^ils  unsteadiness  is  well  documented  in  the  literature.  Stoomel  [107] , 
for  example,  investigated  the  surface  currents  off  Bermuda,  using  drifting 
wireless-telemetering  buoys.  He  found  that  there  was  considerable  irreg¬ 
ularity  in  current  speed  and  direction,  even  during  days  when  the  wind  was 
fairly  steady.  He  calls  this  irregular  motion,  after  Ekman,  a  kind  of  "macro- 
turbulence."  Knauss  [56]  ,  using  neutrally  buoyant  floats,  also  found  "erratic" 
motion  in  deeper  water  below  the  thermocline  (300  to  2l4o  m, ),  Reid  [94]  noted 
that  short-period  variations  In  the  strength  of  the  California  Current  have 
been  observed  of  the  same  order  as  the  strength  of  the  flow  (5-15  cm./sec.). 

The  reasons  tor  this  variability  have  been  sought  in  the  condition  that  there 
are  several  other  forces  besides  wind  and  tides  which,  independent  of  each 
other,  affect  the  surface  currents  simultaneously.  These  are;  (a)  the  "per¬ 
manent  flow"  (either  thermohallne  or  wind-driven  large-scale  circulation); 

(b)  the  changes  of  atmospheric  pressure;  (n)  mass  transport  by  waves;  (d) 
internal  waves;  (e)  the  "modifying  forces"  (corlolls  force,  configuration 
of  the  bottom  and  coast,  permanent  and  inertia  currents);  and  (f)  the  periodic 
tidal  currents.  In  addition,  the  conditions  in  neighboring  areas  (piling  up, 
differences  gradients  in  driving  forces,  etc.)  affect  surface  currents  in  a 
given  area  and  cause  various  eddies  (both  macro  and  meso  scale)  as  well  as  con¬ 
vergences  and  divergences. 

Various  approaches  have  been  used  in  the  past  for  prognostication  of  sur¬ 
face  currents,  but  actual  synoptic  analyses  and  forecasts  have  been  very  lim¬ 
ited  (with  the  exception  of  tidal  currents). 

In  this  manual  an  attempt  is  made  to  present  thuse  facts  sifted  from  the 
literature  and  from  personal  observation  which  have  application  in  analysis 
and  forecasting  of  surface  currents.  The  various  components  of  these  currents 
are  evaluated  separately.  Suggestions  are  given  for  using  these  facts  under 
field  conditions,  and  results  of  large  scale,  synoptic  current  calculations  by 
electronic  computer  are  shown.  In  short,  this  manual  does  not  pretend  to  be  a 
complete  summary  on  ocean  currents  as  such  but  Is  primarily  aimed  at  opera¬ 
tional  current  analysis  and  prediction  on  a  synoptic  time  scale. 

Climatological  surface  current  charts  are  widely  used  for  the  estimation 
of  current  drift.  These  charts,  constructed  from  data  on  ships'  drifts,  are 
available  for  nearly  all  ocean  areas.  Their  accuracy  and  usefulness  have. 
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however,  severe  limitations.  The  data  presented  in  current  atlases  are  useful 
as  a  rough  estimate  of  the  actual  current  at  a  given  time  and  place.  The  re¬ 
ports  upon  which  these  atlases  are  based  were  obteined  from  ships’  logbooks 
end  ere  nothing  more  than  the  difference  between  dead-reckoning  fixes  (from 
ship's  speed  and  direction)  and  navigational  fixes  (astronomical,  Loren,  visu¬ 
al,  etc.). 

Hela  [38]  found  that  if  the  current  displacement  determined  by  astronomical 
fix  is  less  than  2  nautical  miles  per  24  hours,  the  value  in  an  atlas  is  en¬ 
tirely  unreliable.  In  most  cases,  little  weight  can  be  given  to  displacements 
of  5  nautical  miles  or  less  per  24  hours.  James  [51]  came  to  an  analogous 
conclusion)  namely,  that  the  surface  current  charts  can  be  considered  reliable 
and  accurate  as  mean  values  if  there  are  at  least  5.  observations  from  3,600 
square  nautical  miles  and  the  drift  is  more  than  5  miles  per  24  hours.  There 
is  obviously  a  real  need  for  synoptic  current  charts,  particularly  in  offshore 
areas  (the  data  In  some  coastal  waters  are  more  reliable). 

Numerous  detailed  current  measurements  taken  with  meters  of  various  types 
are  now  available  for  limited  areas.  These  have  been  useful  in  building  up 
our  knowledge  and  theories  on  surface  currents.  Most  cf  these  measurements 
have  been  made  relatively  near  the  coast,  and  coastal  and  tidal  effects  limit 
their  use  for  drawing  conclusions  on  offshore  currents.  Many  of  the  current 
measurement  series  are  also  too  short  for  proper  evaluation.  Measurements  of 
actual  surface  currents  offshore  are  extremely  difficult  to  make.  Detailed 
discussion  on  current  measurement  and  various  meters  used  can  be  found  in  ex¬ 
cellent  summaries  by  Wiegel  and  Johnson  [125]  and  Paquette  [88). 

Several  theoretical  discussions  of  surface  currents  in  the  oceans  are 
also  available.  Most  of  the  theoretical  models  are  concerned  with  large  scale 
circulation  and  fall,  therefor*;  outside  the  scope  of  the  present  manual.  The 
reader  is  Termed  to  surcmariea  on  large  scale  circulation  in  the  ocean  by 
Stomnel  [108J ,  Bowden  [10]  and  Robinson  [95]  .  Theoretical  background  and  some 
descriptive  material  on  currents  can  be  found  in  textbooks  such  as  Defant  [20] , 
Sverdrup,  Johnson  and  Fleming  [111] ,  Proudman  [91] ,  Stonmel [109]  and  the  "Bibli¬ 
ography  on  Waves,  Currents  and  Swell"  by  the  American  Meterological  Society  [l] . 

It  has  beeii  customary  to  treat  currents  according  to  either  Lagrange’s 
approach  (following  the  path  of  a  given  fluid  particle)  or  Euler’s  approach 
(describing  flow  characteristics  at  a  given  point).  Neither  of  these  is  fol¬ 
lowed  consistently  in  this  manual.  In  most  cases,  current  is  considered  as  the 
flow  of  water  with  a  defined  average  speed  and  direction  In  a  given  time  inter¬ 
val  (synoptic  period).  However,  a  distinction  is  sometimes  made  between  "mean 
current"  —  which  is  considered  as  a  current  set  with  given  average  speed  and 
direction  in  a  time  interval  greater  than  the  synoptic  period  —  and  "actual 
currenf"  —  which  is  considered  the  current  set  (with  average  direction 
and  speed)  in  a  given  spot  and  time  (time  interval  being  24  hours  or  shorter). 

Although  the  present  stage  of  knowledge  on  ocean  currents  is  not  perfect, 
the  available  material  permits  the  initiation  of  synoptic  analyses  and  fore¬ 
casts.  These  analyses  and  forecasts  are  required  as  inputs  in  other  synoptic 
oceanographic  analyses  and  forecasts  as  well  as  for  direct  application  per  sc. 
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Although  the  variability  of  currents  is  considerable,  attempts  sho'iid  be 
made  to  evaluate  quantitatively  the  causes  of  these  variations.  This  car*  best 
be  done  through  forec-sating  and  daily  verification.  Forecasts  have  little 
value  if  they  cannot  be  verified  by  actual  synoptic  analyses.  Due  to  recent 
developments  in  other  fields  of  oceanographic  analysis  and  prediction,  both  di¬ 
rect  and  indirect  verification  are  now  possible  (see  chapter  5).  Through  these 
verifications  additional  knowledge  is  gained  and  existing  ideas  refined.  It 
is  to  be  expected  that  the  approaches  described  in  this  manual  will  undergo 
considerable  modification  in  the  future. 
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2.  EVALUATION  OF  THE  COMPONENTS  OF  SURFACE  CURRENTS 


The  speed  and  direction  of  the  surface  current  at  a  given  point,  time 
and  depth  below  the  surface  (e.g.,  at  r  depth  of  3  ra.)  \s  can  be  rep*  cscivfceu 
as  the  resultant  of  the  following  components:  ! 
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where: 

Gc  -  permanent  flow  (gradient  current,  and/or  "characteristic  current" ); 

\VC  -  wind  current  (lately  often  erroneously  called  the  "Ekman  current"); 
W,.  -  mass  transport  velocity  by  waves; 

Pi  -  periodic  portion  of  inertia  current; 

P,  -  periodic  portion  of  tidal  current; 

Ac  -  current  caused  by  changes  of  atmospheric  pressure  and  sea  level; 

I  -  the  velocity  and  directional  component,  caused  by  influencing  fac¬ 
tors,  such  as  changing  depth  of  water,  coriolis  force  and  coast  and 
current  boundaries. 

The  direction  and  speed  are  considered  to  be  specified  for  all  components. 
The  present  quantitative  knowledge  of  these  constituent  components  f  applicable 
to  current  predictions,  will  be  evaluated  separately  in  the  following  para¬ 
graphs  . 

2.1  Permanent  Flow 


f 
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Permanent  flow  is  directly  related  to  pressure  gradients  in  the  sea  which 
usually  result  from  density  differences  of  surface  water  columns  between  differ¬ 
ent  locations  ("thermohaline  circulation").  The  permanent  flow  is  maintained 
by  the  great  permanent  and  relatively  constant  wind  systems  ("wind-driven  cir¬ 
culation"  )  which  cause  piling  up  of  water  masses  and  changes  in  sea  level,  re¬ 
sulting  in  pressure  gradients.  It  is  believed  by  some  investigators  that  the 
permanent  flow  is  in  moat  areas  an  "inertia  function"  due  to  prevailing  winds 
and  is  therefore  frequently  referred  to  as  the  "characteristic  current" 

(  Palmer,  [87] ). 

For  any  synoptic  analysis  and  forecast  of  actual  surface  currents  it  is 
necessary  to  separate  the  relatively  persistent,  slowly  changing  (Carruthers 
ll3]  and  Carruthers,  Lawford  and  Veiey  [id] )  permanent  flow  from  other  compo¬ 
nents  of  surface  currents  which  may  vary  hourly  and  daily  and  which  ahould^be 
predicted  using  the  direct,  synoptic  influence  of  their  driving  forces. 

In  subjective  forecasting  (small-area),  the  permanent  flew  should  usually 
be  estimated  using  separation  methods  and  then  be  added  to  other  components 
computed  and/or  estimated  by  the  synoptic  consideration  of  their  driving  forces. 
In  large  scale  numerical  (objective)  forecasts  the  permanent  flow  can  be  com¬ 
puted.;  however,  "separated"  permanent  flow  should  also  be  used  to  verify  th^se 
computations. 

The  basic  method  for  separation  of  permanent  flew  has  been  developed  and 
tested  by  Palmen  [87]  and  Kela  [37] „  They  took  the  currents  associated  with  a 
certain  wind  velocity,  grouped  them  according  to  the  directions  of  the  winds 


1 


and  currents  and  computed  the  resultant  currents  for  each  direction  (see  for¬ 
mulas  for  wind  current  in  section  2,2),  The  permanent  flow  was  obtained  by 
subtracting  the  computed  current  from  the  measured  current.  This  method  of 
separating  out  thf  permanent  flow  yields  in  most  cases  the  so-called  "charac¬ 
teristic  current,"  thus  indicating  that  the  total  surface  currents  are  mainly 
wind-driven. 


In  areas  where  the  data  on  average  monthly  surface  currents  are  not  avail¬ 
able  or  are  unreliable  because  of  too  few  observations,  it  is  necessary  to  ap¬ 
proximate  the  permanent  flow  using  other  information  such  as  sea  level  data, 
hydrographic  (oceanographic)  station  data  and  observed  winds. 


Formulas  have  been  derived  for  computation  of  currents  from  a  pair  of  hy¬ 
drographic  stations  (utilizing  tne  difference  in  the  dynamic  heights  at  the 
stations)  and  for  computation  of  gradient  flow  in  a  <J,  surface. 


The  formula  of  Sandstrom  and  Helland-Hansen  for  computation  of  a  friction- 
less,  stationary,  relative  current  component  cross  the  section  between  two 
stations  is* 


w*  -  w5 
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—11 -  (AD.  -  AD„) 

2  £  co  sin  <f>  A  B 


(2) 


(symbols,  see  appendix).  The  computation  of  AD  (the  dynamic  height  anomaly) 
from  hydrographic  data  is  described  in  H.  0.  Pub.  607  (U.S.  Navy,  Hydrographic 
Office  [117]).  The  relative  current  is  positive  (directed  away  from  the  read¬ 
er)  if  station  A  lies  to  the  right  of  B  and  if  ad.  Is  greater  than  ADQ  . 

A  B 


Equation  (2)  cannot  be  used  in  its  basic  form  for  computations  of  gradient 
currents  in  shallow  water.  However,  Groen  [32]  has  described  an  approximate 
method  for  estimating  slopes  and  currents  in  shallow  water.  lie  replaced  the 
solid  earth  by  a  fictitious  section  of  water  in  which  the  density  distribution 
was  determined  by  putting  the  inclination  of  the  ■„  lines  along  every  horizontal 
(or  isobarlc)  line  equal  to  the  value  projected  horizontally  from  the  bottom 
line.  The  lines  were  then  extended  from  the  actual  water  into  the  fictitious 
vetfrr  by  reference  to  these  inclinations. 


Although  the  original  hydrodynamic  approach,  equation  (2),  was  developed 
for  use  between  two  quasi -synoptic  hydrographic  stations,  its  use  has  been 
widened  (but  its  limitations  often  neglected). 


Mapping  of  dynamic  height  topography  in  relation  to  a  reference  level  (e.g., 
1000  m. )  is  assumed  to  give  a  general  picture  of  circulation  if  enough  quasi - 
synoptic  hydrographic  stations  are  available  from  the  area  (see  fig.  2.1)  Ap¬ 
proximate  current  speeds  are  obtained  by  measuring  the  distances  between  con¬ 
tours. 


Montgomery  [73]  concluded  that  a  logical  method  of  deducing  oceanic  flow 
patterns  (permanent  flow)  from  temperature,  salinity  a. id  oxygen  observations 
is  to  chart  these  properties  for  a  surface  of  constant  potential  density.  His 
£"P“cs5ion  for  the  gradient  flew  in  a  ,  surface  is: 


37  AHg 

sin  4>  £ 


(3> 
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where  Hg  (the  geostrophic  potential)  is: 

<t>n  ~  *0P* 

The  following  Just  criticisms  have  bean  raised  concerning  the  application 
of  equations  (2)  and  (3): 

(a)  The  requirement  that  the  stations  in  the  section  are  quasi-synoptic 
is  seldom  fulfilled. 

(b)  The  computation  of  AD  requires  an  assumption  of  a  "level  or  layer 
of  no  motion,"  which  might  exist  but,  at  best,  is  difficult  to  assess  correct¬ 
ly  (severed,  theoretical  approaches  for  estimation  of  this  layer  exist,  e.g., 
Befant[20]).  This  layer  is  generally  assumed  to  be  between  500  and  2000 
meters.  The  results  obtained  from  equation  (2)  also  depend  upon  the  spacing  of 
stations . 

(c)  The  equation  cannot  be  used  near  the  equator.  (A  review  of  the  calcula¬ 
tion  of  ocean  currents  at  the  equator  has  been  made  by  Arthur  [3] ) . 

Sane  of  the  gradient  currents,  such  as  the  equation  counter-current  in 
the  central  Pacific,  are  remarkably  predictable  in  latitude  and  general  physi¬ 
cal  characteristics  (Austin,  Stroup  and  Rinkel  [41  ).  However,  in  most  areas 
dynamical  calculations  with  equations  (2)  give  only  a  rough  idea  of  the  cur¬ 
rents.  Saelen  [98]  has  made  extensive  comparison  of  geostrophlc,  computed  and 
measured  currents  in  the  Norwegian  Sea.  Often  he  found  only  a  qualitative 
agreement.  Saelen  (op.  cit. )  furthermore  found  that  the  conception  of  the 
circulation  of  surface  and  deep  waters  in  the  Norwegian  Sea  as  a  slow  and  regular 
drift  does  not  hold,  and  currents  at  all  levels  vary  rapidly  over  short  periods 
of  time. 

As  hydrographic  stations  are  seldom  synoptic,  tidal  fluctuations  may  lead 
to  considerable  errors  in  distribution  of  dynamic  heights,  Defant  [19]  has  of¬ 
fered  a  method  for  elimination  Of  "tidal"  errors  by  preparation  of  dynamic 
charts  for  basins  with  simple  tidal  circulations,  and  Zaitsev  [133]  suggested 
a  modification  of  this  method  for  seas  with  complex  tidal  circulation,  utiliz- 
ing  the  co-tidal  and  co-range  lines  with  coastal  tidal  data.  For  these  cor¬ 
rections  the  following  data  are  needed:  (a)  hydrographic  station  data;  (b)  co- 
tidal  chart  (usually  theoretically  computed  M2  component)  of  the  area  where 
co-range  lines  are  indicated;  (c)  tidal  curve  for  the  period  of  the  hydLrograhic 
stations  from  one  of  the  closest  unsheltered  representative  coastal  tide-gauge 
stations  (reference  station).  The  suggested  correction  procedure  for  an  indi¬ 
vidual  station  is  —  "The  actual  station  time  is  noted  on  the  tidal  amplitude 
curve  or  the  reference  station.  The  co-tidal  difference  between  the  hydrograph¬ 
ic  station  and  the  reference  station  is  estimated  from  the  co-tidal  chart,  and 
the  time  difference  is  transferred  to  the  amplitude  curve  of  the  reference  sta¬ 
tion  (either  to  the  right  or  left  of  the  noted  actual  station  time).  From  this 
point  the  tidal  amplitude  is  read  from  the  reference  station  curve.  This  ampli¬ 
tude  is  corrected  with  the  ratio  of  amplitude  differences  between  reference  sta¬ 
tion  and  hydrographic  station,  and  the  resultant  amplitude  is  used  to  correct 
the  actual  sampling  depths  — ."  This  method  does  not  account  for  internal  tidal 
effects,  which  are  thought  to  be  considerable  (FJeldatad  [26])  but  have  not  yet 
been  evaluated  quantitatively. 
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Worthington  [129]  found  that  dynamic  computations  of  the  Gulf  Stream,  using 
closely  spaced  stations.,  indicate  that  the  cross-current  slope  of  the  iaobarie 
surfaces  is  uneven,  resulting  theoretically  in  three  or  four  aonee  of  high  ve¬ 
locity  separated  by  za ne?  of  relatively  lower  velocity.  He  suggested  that  the 
unevenness  of  the  slopes  of  the  itobaric  surface  is  related  to  the  surface  tem¬ 
perature  discontinuities  observed  in  the  Gulf  Stream  by  von  Arx  fix d  Richardson 
(1953) >  end  that  the  stream  is  composed  of  overlapping,  discontinuous  currents 
at  all  levels,  If  this  suggestion  is  correct,  it  would  be  virtmliy  impossible 
to  determine  a  "level  of  no  motion"  and  apply  hydrodynamical  methods  to  estima¬ 
tion  of  currents  in  this  region,  (see  also  fig.  g.2). 

In  view  of  the  limitations  described  so  far  In  this  chapter,  it  becoaee 
obvious  that  the  conventional  hydrodynamics!  approach  is  unsuitable  for  numer¬ 
ical  synoptic  analyses  eu.d  forecasting  of  the  permanent  flow.  The  main  draw¬ 
backs  are  the  following*. 

(a)  The  dynamic  method  gives  only  s,  qualitative  picture  of  the  currents, 
which  can  vary  considerably  over  short  periods. 

(b)  It  is  nearly  impossible  to  obtain  synoptic  density  distribution:  the 
only  good  approximation  must  come  frets  synoptic  thermal  structure  analyses. 

(c)  The  dynamic  computations  give  the  total  current,  but  only  the  perma¬ 
nent  flow  component  is  desired,  since  other  components  setan  to  be  more  easily 
attacked  by  data  on  their  driving  forces. 

It  is  therefore  suggested  that  the  beet  approach  for  numerical  synoptic 
computation  of  permanent  flew  is  to  compute  it  from  synoptic  thesraal  and  clima¬ 
tological  salinity  structure  (tuned  with  appropriate  constants )  based  on  ser-sra 
tion  data  at  a  number  of  localities.  Such  an  approach  is  in  use  at  the  Fleet 


in  Late  Winter  1958.  Positive  Values  Give  Current  to  E  or  SE  and  Negative  Values 
(gray  areas)  Give  Current  to  W  or  NW.  (Dynamical  Computations;  after  Koslowski, 
[58]  modified) 


Numerical  Weather  Facility  (Hubert  [43]).  Indirect  support  for  this  approach 
can  be  found  in  a  few  specific  studies:  Yeaui [131]  found  that  there  is  a  close 
correlation  between  the  distribution  of  temperature  below  the  thermoeline  and 
dynamic  depth  anomalies.  He  constructed  a  simple  empirical  relationship  by 
plotting  temperatures  versus  computed  dynamic  depth  anomalies,  thereby  shorten¬ 
ing  the  dynamic  computations.  This  relation,  once  established  for  a  given  lo¬ 
cation,  has  been  assumed  to  be  valid  for  a  wider  area. 

Neglecting  salinity,  one  can  apply  the  well-known  meteorological  thermal 
wind  relationship  in  the  ocean  if  one  knows  the  mean  temperature  of  the  layer 
between  the  surface  and  some  deeper  level  of  alow  current  velocity  (?.eiected 
as  a  level  with  negligible  annual  temperature  change).  The  characteristic  sur¬ 
face  current  is  then  given  by: 

W  =  -III-  7Tx|K  (4) 

f  T 

where  g  is  gravitational  acceleration,  f  la  the  coriolls  parameter,  T  is  the 
mean  temperature  above  the  level  of  assumed  zero  current.  As  l.a  the  depth  to 
zero  current  and  |K  is  the  unit  vertical  vector. 

The  determination  of  a  representative  mean  temperature  (T)  is,  of  course, 
a  critical  factor  in  this  problem.  The  temperature  structure  of  the  ocean  is 
certainly  not  constant,  particularly  closer  to  the  surface;  therefore,  synoptic 
temperature  fields  should  be  used  if  possible.  The  only  place  where  sufficient 
data  are  available  for  reliable  analysis  on  a  daily  hauls  is  at  the  surface. 
(Fleet  Numerical  Weather  Facility  Sea-Surface  Temperature [SST] analyses . )  In 
order  to  Include  a  part  of  the  deeper  temperature  structure,  the  SST  field  is 
combined  with  a  climatological  field  at  200  meters  depth  to  obtain: 

T  ”  TaJc  +  K£  T200  (5) 


Finally,  this  field  cen  be  modified  empirically  In  areas  where  salinity 
considerations  are  known  to  be  important  (Oyashlo,  Greenland,  Labrador  cur¬ 
rents).  This  in  effect  corrects  the  ocean  temperatures  for  salinity  much  as 
the  meteorologist  corrects  atmospheric  temperatures  for  moisture  content  when 
he  uses  the  concept  of  "virtual  temperature."  The  Justification  of  the  use  of 
200  m.  temperatures  and  details  of  the  nunerical  methods  are  given  in  chapter  4, 

Some  additional  hydrodynamic  approaches  for  estimation  of  currents  from 
sea  level  and/or  interface  slope  are  given  in  sections  2.3  and  2.4.  These 
approaches  are  of  use  where  synoptic  sea  level  records  are  available. 

2.2  Wind  Currents  and  Mass  Transport  by  Waves 

The  wind  is  one  of  the  main  driving  forces  of  surf Me  currents,  and  many 
current  systems  are  mostly  wind-driven.  Wind  is  directly  or  indirectly  one  of 
the  main  causes  of  temperature  fluctuations  in  the  sea  and  thus  also  affects 
the  permanent  flow. 

In  numerous  cases  a  simple  linear  relation  between  the  wind  and  current 
velocities  has  been  found  and/or  used  In  the  past.  The  current  speed  has  been 
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found  to  vary  from  1.2%  to  4%  of  surface  wind  speed.  'rh!i=  lb  evidence  that  the 
factor  decreases  with  increasing  wind  speed,  an  observation  which  has  caused 
some  confusion  and  led  to  much  discussion  concerning  a  possible  "critical" 
wind  velocity.  There  is  also  some  evidence  that  fetch  length  and  wind  duration 
have  an  effect  on  wind-currents;  however,  these  factors  have  not  been  properly 
evaluated. 

A  simple  relation  can  be  used  to  approximate  wind  drift  from  climatologi¬ 
cal  data  (monthly  resultant  wind  speeds): 

W  =  KjV  (6) 

The  factor  K,  varies  between  0.014  and  0.025.  The  value  of  0.02  has  been 
found  most  appropriate  for  a  rough  estimation  of  the  sum  of  average  wind  and 
"wave"  current  in  the  surface  layer  (at  about  3  ro.  depth).  Although  Thorade 
(from  Defant  [20]  )  assumed  K,  to  be  a  function  of  latitude,  there  is  no  solid 
theoretical  reason  or  empirical  evidence  for  this. 

The  formula  of  Witting [127]  ,  based  on  more  than  5000  observations,  has  been 
generally  found  to  present  best  the  relation  between  the  actual  wind  and  act¬ 
ual  current  velocities; 

w  =  K^fyT  (7) 

where  W  is  current  velocity  (cm. /sec.)  and  V  is  wind  velocity  (m./sec.).  The 
factor  K2  varies  slightly  with  depth  (see  section  3*1) •  As  an  average,  the 
value  of  3.8  can  be  used  if  the  mass  transport  by  waves  ("wave  current")  is 
estimated  separately  and  an  average  current  speed  of  the  mixed  layer  is  re¬ 
quired.  (If  mass  transport  by  waves  is  included  in  this  formula,  the  factor 
is  between  4  and  5). 

According  to  Ekman's [25]  theory,  the  direction  of  the  wind  current  at  the 
surface  is  45*  to  the  right  in  the  northern  hemisphere,  and  the  current  speed 
decreases  vith  depth  while  the  deflection  becomes  larger  until,  at  a  "depth  of 
frictional  Influence,"  the  direction  of  the  current  is  opposite  to  that  at  the 
surface.  Although  Bkman* s  theory  has  been  useful  in  many  Instances,  it  has  not 
been  verified  in  detail  in  the  field.  Ekman's  theory  cannot  be  applied  to 
practical  situations,  because  the  conditions  under  which  the  problem  was  solved 
do  not  exist  in  nature  (Belinsky  and  Glagoleva, [8]  ). 

Recent  investigations  indicate  that  the  wind  current  in  offshore  areas  in 
middle  latitudes  deflects  in  general  less  than  20*  to  the  right  of  the  wind. 
(Bowden,  [1953]  -  18*;  Gaul  [30]  -  15*;  Lisitzin  [68]  -  12*;  Stomnel  [107]  -20*; 
see  also  Manner  [72]  ).  Furthermore/,  it  has  been  found  that  the  current  de¬ 
flection  depends  also  on  the  wind  speed  (Manner  [72]  ) ,  Therefore  the  empirical 
formula  of  Hela  [37],  derived  for  conditions  in  the  Baltic  Sea,  merits  further 
testing  in  different  latitudes  and  in  the  open  ocean; 

8  -  34°-7.5°yV  /Al 


where  y9  is  the  deflection  in  degrees  to  the  right  of  the  wind,  and  V  ia  the 
climatological  mean  wind  speed  (m./sec.). 
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lt  the  wind  Is  weak,  other  current  components  can  prevail,  and  the  current 
movement  is  "erratic"  in  relation  to  the  wind.  At  lower  wind  velocities  the 
angle  of  deflection  is  somewhat  larger,  shout  25*  to  30*  in  higher  latitudes, 
but  converges  rapidly  between  wind  speeds  of  4  to  6  knots  to  about  20* 
(Shulekin), 


Tentatively,  the  following  deflection  angles  are  suggested  for  winds  of 
6  knots: 

Latitude  Deflection  (degrees) 
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With  lower  winds  the  deflection  angle  should  be  increased  by  4*  to  12*. 

Some  observations  of  the  behavior  of  surface,  currents  in  oceanic  areas 
indicate  that  the  wind  current  is  not  always  determined  by  the  wind  prevail¬ 
ing  over  a  given  location  but  rather  by  moving  wind  fields  (and  especially 
by  stronger  winds)  over  a  larger  area  creating  a  relatively  complex  picture. 


The  relation  between  sea  level  changes  (caused  by  wind  set-up)  and  cur¬ 
rent  has  been  rather  extensively  studied  by  Palmen  and  Lisitzin.  Their  re¬ 
sults  might  bo  of  use  for  estimation  of  currents.  A  formula,  based  on 
Torricelli's  theorem  and  extensively  tested  by  Lisitzin  (e.g., [69]  ),  is: 


(9) 


The  static  effects  of  the  atmospheric  pressure  differences  must  be  eliminated 
beforehand]  the  slope  of  the  surface  (AHW)  is  given  in  cm.  per  100  km., wind 
(V)  in  m./see.  and  depth  of  water  (d)  in  m.w 


Mass  transport  by  waves  has  been  quantitatively  accounted  for  in  very  few 
earlier  surface  current  hindcasting  or  forecasting  attempts  but  has  been 
indirectly  (and  even  unconsciously)  included  in  empirical  wind  current  formu¬ 
las.  This  inclusion  is  advisable  also  in  the  future  because  existing  uncer¬ 
tainties  in  quantitative  evaluations  do  not  Justify  a  separate  computation. 

In  certain  cases,  especially  in  nearshore  problems  (e.g.,  longshore  currents 
and  currents  in  reef  channels,  etc.),  it  is  necessary  to  account  separately 
for  mass  transport  by  waves  and  for  particle  velocities  in  waves  at  different 
depths  below  the  surface. 


Sane  laboratory  experiments  on  the  mass  transport  by  waves  and  their 
theoretical  treatment  are  available  in  the  literature  (e.g.,  Masch [73l ; 
Longuet-Higgina  (?ol ) .  Unfortunately,  laboratory  experiments  in  wave  tanks 
(for  the  investigation  of  mass  transport)  do  not  reproduce  conditions  found 
in  nature.  The  transport  is  quite  different  in  the  beginning  of  the  experi¬ 
ment  when  the  wave  generator  is  started  (forward  transport  on  the  surface 
only).  After  a  relatively  short  time,  a  circulation,  determined  by  the  tank 
boundaries,  is  established  in  the  wave  tank  (backward  transport  on  the  surface 
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In  the  case  of  a  mechanical  wave  generator,  vice  versa  by  "blower"  generator, 
forward  transport  at  same  Intermediate  depth  or  on  the  bottom,  depending  on 
the  relations  between  the  wave  length  and  the  depth  of  the  water). 


According  to  theory,  Qerstner  waves  are  trocholdal  and  without  mass  trans¬ 
port  velocity;  Stokes  waves,  however,  have  a  transport  velocity  {lamb  [63]; 
Pershin  and  Voznessensky  [89] ;  Sverdrup  and  Munk  [112]).  A  theoretical  formula 
for  mass  transport  velocity  is  usually  given  in  the  following  forms 


W 


m 


(10) 


This  formula  can  be  used  with  certain  assumptions: 


Lj/3  -  50-^H j^j  ;  C  «  1.25  L1/3  (ll) 

Masch  [73]  found  that  equation  ( 10)  gave  values  too  small  for  mass  transport 
as  compared  to  the  results  of  laboratory  experiments.  He  concluded  that  the 
following  equation  ( 12)  described  the  mass  transport 


W  -  4n2  (  Hw  )*  c  <12) 

m  \  '  o  i 

Li/3  Li/a 

The  mass  transport  by  waves  decreases  rapidly  with  depth,  and  the  effect 
on  the  drift  of  a  ship  depends  on  the  ship's  draft.  There  is,  unfortunately, 
no  verification  of  the  theoretical  treatment  and  laboratory  measurements 
-trough  observations  at  sea.  In  wave  tanks,  equation  (12)  accounted  for  both 
wind  current  and  mass  transport  by  waves.  Further  tests  of  this  equation  in 
the  sea  are  necessary  (see  fig.  2.3). 


MASS  TRANSPORT  VELOCITY,  U  (ft. /sec.) 

Figure  2.3  Velocity  of  Mass  Transport  at  Various  Depths  as  a  Function  of  Wave  Period 
(from  Beach  Erosion  Board,  1942) 

As  the  direction  of  mass  transport  is  in  most  cases  approximately  the 
same  as  the  wind  current  component,  and  as  mass  transport  is  of  importance  only 
close  to  the  surface,  practical  surface  current  analyses  and  forecasts  can  ac- 
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count  for  this  term  by  an  empirical  increase  of  the  total  effect, 

2.3  Tidal  and  “Hydraulic”  Currents 

A  knowledge  of  tidal  currents  is  of  great  importance  for  navigation  in 
coastal,  waters,  especially  in  inlets.  Therefore,  predictions  of  tidal  currents 
are  made  for  a  number  of  important  locations.  The  interpolation  between  the 
prediction  localities  is,  however,  in  many  cases  rather  uncertain.  Tidal  cur¬ 
rent  analysis  and  prediction  for  oif shore  waters  have  been  greatly  neglected, 
and  it  has  been  generally  assumed  that  tidal  effects  are  weak  over  deep  water. 
Latest  measurements,  however,  indicate  that  tidal  currents  can  have  consider¬ 
able  magnitude  even  in  mid-ocean .  For  a.  24-hour  prediction  the  tidal  currents 
might  be  neglected  in  most  cases.  However,  in  detailed  interpolation  of  these 
forecasts,  the  tidal  currents  must  be  accounted  for.  Some  specific  aspects  of 
tidal  currents  are  pointed  out  below: 

(a)  The  tidal  currents  are  usually  rotary  offshore  and  reversing  in  coast¬ 
al  waters  (see  Hebard  [36]).  The  tidal  ellipses  get  narrower  when  one  approach¬ 
es  the  coast. 

(b)  In  the  northern  hemisphere  the  sense  of  rotation  is  usually,  but  not 
always,  dextral  (clockwise). 

(c)  In  general,  the  flooding  current  runs  in  the  direction  of  the  movement 
of  the  tidal  wave  (the  direction  can  be  ascertained  from  a  co-tidal  line  chart). 

(d)  The  rising  tide  is  not  everywhere  synonymous  with  the  flooding  current, 
and  falling  tide  does  not  necessarily  mean  ebbing  current. 

(e)  High  tide  does  not  necessarily  mean  strong  current,  and  a  strong  cur¬ 
rent  may  accompany  a  small  tide.  It  is,  therefore,  always  necessary  to  dis¬ 
tinguish  clearly  between  tide  and  tidal  currents.  The  relations  between  the 
two  are  not  simple  and  not  the  same  in  all  locations. 

(f)  The  behavior-  of  tidal  currents  of  otherwise  similar  tides  is  affected 
by  the  difference  in  location  of  amphidromir.  points  of  different  constituents. 

(g)  Semidiurnal  tide -producing  forces  cause  proportional  changes  in  tides 
and  currents,  but  daily  forces  affect  tho  current  only  half  as  much  as  the  tides. 

(h)  The  ebb  usually  lasts  longe1"  than  the  flood  in  coastal  waters.  How 
this  can  apply  to  rotary  or  nearly  rotary  currents  offshore  is  ;rot  known. 

(i)  Some  recent  recordings  of  tidal  currents  in  offshore  waters  indicate 
that  the  14-day  periodicity  is  a  very  prominent  feature.  (Mossby,  personal 
communication. )  This  14-day  periodicity  has  been  pointed  out  by  Laevaatu  [60] 
in  the  fluctuations  of  themocline  depths  as  well  as  in  fluctuations  of  sea 
level . 

( J )  The  topography  of  the  coasted,  area  influences  the  tidal  currents  con¬ 
siderably  (see,  e.g.,  Neumann  [83])  causing  eddy  formation  behind  headlands  and 
in  bays  (Laevastu,  Avery  and  Cox  [61]). 

(k)  The  tidal  current  ellipses  are  narrow  and  the  currents  strong  in  chan¬ 
nels  and  on  submarine  slopes. 

(l)  The  tidal  current  turns  sooner  in  shallow  water  (close  to  the  coast) 
than  farther  offshore. 

An  example  of  the  behavior  of  tidal  currents  in  relation  to  the  tidal  ele¬ 
vation  curve  in  the  Hawaiian  Islands  area  Is  shown  on  figure  2.4. 


In  view  of  the  description  of  tidal,  currents  above,  it  can  be  concluded 


Recordings  of  Tidal  Currents  in  Two  Locations 
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that  the  only  sufficiently  accurate  method  for  prediction  of  tidal  currents  in 
coastal  waters  Is  to  record  the  tidal  currents  at  a  number  of  points  along  the 
coast  for  at  leut  25  hours  (preferably  15  to  29  days),  to  analyze  the  changes 
of  directions  and  speeds  in  relation  to  tidal  elevation  curves  and  then  make 
predictions  on  the  basis  of  these  analyses,  (See  relation  between  tidal  ampli¬ 
tude  and  currents,  Marser  [72],)  The  currents  between  points  of  observations 
must  be  interpolated.  The  accuracy  of  this  interpolation  usually  increases 
with  decreasing  distance  between  the  points.  If  necessary  and  possible,  other 
current  components  should  be  subtracted  during  the  analyses  and  not  included 
in  the  prediction  of  the  tidal  components  themselves.  Such  predictions  (based 
on  empirical  knowledge  in  a  given  location  of  the  tidal  currents  in  relation 
to  tides)  remain  valid  indefinitely. 

At  times  harmonic  analysis  and  prediction  of  currents  can  also  be  made. 
This  is  usually  a  standard  procedure  for  the  preparation  of  tidal  current 
tables.  The  methods  of  harmonic  analyses  and  prediction  of  tides  and  tidal 
currents  are  described  in  manuals  such  as  U.S.  Coast  and  Geodetic  Survey 
(Schureman  [1011),  "Manual  of  Harmonic  Analyses  and  Prediction  of  Tides"  or  in 
Admiralty  of  U.K.  (Suthons  [110])  "Admiralty  Tidal  Handbook." 

In  coastal  waters  with  complex  coastal  topography  and  in  narrows,  many 
so-called  "hydraulic  effects"  (differences  in  sea  level)  affect  the  tided  cur¬ 
rents.  The  "hydraulic  currents"  in  narrows  are  usually  proportional  to  level 
differences  and  can  be  computed  and/or  predicted  if  the  level  differences  at 
two  locations  and  their  changes  with  time  are  known.  The  "hydraulic  currents" 
are  known  to  exist  in  large  areas  in  large  archipelagos.  If  detailed  current 
predictions  in  such  areas  are  needed,  these  "hydraulic  currents"  can  be  ac¬ 
counted  for  with  hydrodynamic  computations,  using  sea  level  records.  In  field 
conditions  and  for  occasional  predictions,  short-cut  methods  (utilizing  sea 
level  slope)  have  been  used  with  success  in  the  past. 

If  the  atmospheric  pressure  were  uniform  and  there  were  no  "stationary" 
piling  up  of  waters  by  wind  action  along  the  coast,  the  sea  surface  could  be 
considered  as  an  lsobarlc  surface,  and  the  gradient  equation  might  be  used  to 
connect  the  inclination  of  the  sea  surface  (tan  i)  tc  the  surface  current: 

ta„  i  .  2?~gi-n*  W  (13) 

g 

In  the  northern  hemisphere,  in  the  direction  of  the  current,  the  lsobarlc  sur¬ 
face  would  then  slope  upward  to  the  right. 

As  shown  by  Sands trom  (1903),  (from  Hela[39]),  there  must  exist,  roughly, 
the  following  relationship  between  the  transversal  slope  of  the  sea  surface  (i) 
and  the  mean  current  speed  (W)  : 

i  -  2  <ii  sin  <fi  w  (lJ*) 


In  this  formula,  AH,  is  the  transversal  height  difference,  and  Be  is  the 
breadth  of  the  current.  Hela  [37]  found  the  following  empirical  relation  be¬ 
tween  the  water  height  difference  in  cross  currents  and  current  speed  in  the 
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Strait  of  Florida: 


AH,  =  0.0152W2 


where  AM  is  in  cm.  and  >¥  ia  nautical  miles  per  day. 

According  to  Torricelli's  theorem  the  current  speed  can  also  be  related  to 
the  longitudinal  height  difference  (Heia  [37]  ): 

Wj  =  \/'2jAH'£  (1 C) 

Mar gules'  equation  refers  to  the  inclination  of  the  interface  (pycnocline) 
in  relation  to  currents ; 

2  m  sin  (p i W,  -  psff. )  .  . 

tan  y  - -  1 J- 1 1 

&(p2  -  Pi) 

In  a  steady  equilibrium  state,  the  slope  of  the  internal  boundary'  surface 
in  a  two-layered  sea  will  be  greater  than  that  of  she  physical  sea  level  in 
the  ratio: 

Pi}  (*2- Pi)  (18) 

The  interface  slope  (i)  should  then  be  given  by: 

i  „ - i -  (19) 

’ .  ...  '^pi  ~  pi)  hi 

Dietrich  [22]  has  adopted  the  gradient  current  approach  and  the  sloping 
interface  relation  for  computation  of  currents  in  a  lower  layer  along  a  sloping 
bottom  with  the  following  formula: 

w  _  ./'**** s (P*  ~  Pi)  _  (pr>\ 


The  problems  of  prediction  of  tidal  currents  in  offshore  waters  should  be 
susceptible  to  numerical  attack  utilizing  strictly  hydrodynamic  approaches. 
However,  in  the  absence  of  such  predictions,  a  subjective  "field"  procedure  is 
suggested  below  (this  approach  requires  verification). 

If  required,  the  tidal  velocity  component  of  a  current  in  offshore  waters 
could  be  approximately  accounted  for  in  short  tern  forecaste  (for  example,  for 
trawling  purposes)  using  the  tidal  manuals  for  coastal  areas  and  the  co-tidal 
and  co-range  charts.  The  suggested  simple,  very  tentative  procedure  for  estima¬ 
tion  of  tidal  currents  in  offshore  areas  is  as  follows: 

(a)  Assume  that  the  tidal  ws,ve  in  offshore  area  is  a  progressive  wave, 
which  means  that  maximum  current  velocities  occur  at  the  creSt  and  trough  of 
the  tidal  amplitude  curve  for  a  given  location,  and  about  half  of  this  velocity 
results  between  these  points  (i.e.,  a  1:2  tidal  ellipse). 

(b)  Ascertain  the  stage  of  the  tide  and  its  amplitude  in  a  given  location 
and  time  from  the  co-tidal  chart  of  Ms  tide  and  closest  coastal  reference  sta- 
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tior  along  an  open  unprotected  coast. 

(c)  Assume  that  the  current  is  a  function  of  depth 

\/h7 

where  W{  is  the  tidal  current  at  the  reference  station,  and  h,  is  the  depth 
at  any  location.  Multiply  this  result  with  the  factor  of 

\/A 7 

/a; 

where  A,  arid  A2  are  the  tided  ranges  at  the  location  in  question  and  at  the 
reference  station  respectively.  (Assumption  (c)  is  a  very  tentative*  one.) 

. (d)  Assume  that  at  the  tidal  wave  crest  the  current  runs  in  the  direc¬ 
tion  of  ths  movement  (progress)  of  the  tidal  wave  (direction  from  co-tidal 
chart),  and  at  the  trough  it  is  opposite  to  the  movement  of  the  tidal  wave 
(i.e.,  into  the  direction  from  which  the  wave  came).  Assume  that  the  rota¬ 
tion  of  the  current  is  clockwise. 

It  should  be  underlined  that  the  above  suggested  procedure  is  largely  hy¬ 
pothetical  and  not  sufficiently  tested. 

Earlier  attempts  with  different  theoretical  procedures  indicate  that  the 
computed  and  observed  velocities  of  tidal  currents  do  not  agree  well  in  the 
open  ocean.  However,  for  a  24-hour  prognostic  current  forecast  for  navigation, 
the  tidal  currents  could  be  ignored  because  the  resultant  tidal  transport  for 
this  period  will  be  ciose  to  0. 

2.4  Inertia  Currents  and  Other  Current  Components 

When  there  is  no  apparent  resultant  force  (except  the  deflecting  force  of 
the  earth's  rotation)  acting  upon  a  quantity  of  water  in  a  given  locnlity,  the 
current  is  referred  to  as  an  inertia  current.  The  primary  force  causing  the 
inertia  is  usually  the  past  prevailing  and/cr  stronger  winds  in  a  given  area, 
sometimes  at  some  distance. 

The  periodic  changes  of  this  inertia  current  follow  the  theory  (period  and 
circle  of  inertia),  although  the  forecasting  of  the  phase  at  a  given  time  and 
locality  is  nearly  impossible  at  present. 

The  period  of  inertia  is: 


rji  _  1  TT  ^ 

p  2  o>  sin  <f> 


(21) 


T  is  12  hours  at  the  poles,  2k  hours  at  30*  latitude  and  infinity  at  the  equator. 


The  radius  of  inertia  is: 


W 


r 


2  oj  sin  d> 


(22) 
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Available  investigations  on  inertia  currents  show  that  they  affect  the 
speed  of  the  permanent  flow  at  the  given  locality  by  periodic  fluctuation. 

At  present,  no  method  exists  for  making  a  hindcast  or  forecast  of  these  fluc¬ 
tuations;  however,  large  and  fast  inertia  eddies  are  often  observed  (e.g,, 
Defsnt[17],  Barkley,  personal  c oianuni cation).  It  might  be  possible  in  the 
future  to  treat  theoretically  the  inertia  cf  a  surface  current  system  (created, 
for  example,  by  a  passing  storm)  and  its  response  to  new  modifying  forces 
(winds)  by  utilizing  known  relations  of  the  moments  of  inertia  and  kinetic 
energy.  It  might  also  be  possible  to  account  for  inertia  currents  (caused  by 
strong  wind  fields)  in  continuous  subjective  analyses  of  large  areas  when  more 
knowledge  is  accumulated  from  synoptic  current  analyses  and  forecasts  and 
their  verification. 


Theoretically,  it  can  be  expected  that  in  the  open  ocean  an  atmospheric 
pressure  change  will  cause  a  wave  in  the  pycnocllne  rather  than  a  true  sur¬ 
face  current.  However,  there  must  be  seme  current  and  mass  transport  connected 
with  this  wave.  No  direct  empirical  evidence  is  available  on  the  current  com¬ 
ponent  caused  by  the  change  of  atmospheric  pressure.  Recently,  however,  Donn 
and  McGuiness  [23]  have  shown  that  under  certain  conditions  there  can  be  a  cou¬ 
pling  of  air  pressure  and  ocean  waves  which  can  produce  sea  level  oscillations 
100  times  greater  than  the  barometric  equilibrium  wave.  There  is  also  the 
possibility  that  these  "atmospheric  pressure  waves"  are  coupled  vith  the  in¬ 
ertia  currents,  especially  when  the  inertia  currents  do  not  show  dissipation 
with  time  (Defant  [17] ). 

Internal  waves,  with  tidal  periods,  may  cause  intermittent  currents  in 
the  surface  layers.  The  currents  could  be  included  in  a  prediction  of  tidal 
currents  when  accurate  tidal  current  prediction  techniques  are  available  for 
offshore  areas. 


Some  simple  theoretical  relations  for  the  period  and  velocity  of  the  prop¬ 
agation  of  internal  waves  are  available  for  use  in  checking  the  importance  of 
these  waves  in  respect  to  currents: 

(23) 


(24) 


At  present,  no  tested  method  exists  for  predicting  current  component 
caused  by  the  change  of  atmospheric  pressure  or  by  internal  waves. 

The  current  component  caused  by  coriolis  force  does  not  appear  in  the  for¬ 
mulas  given  earlier.  This  is,  however,  not  strictly  necessary,  because  the  ra¬ 
tional  formulas  are  approximations  and  are  usually  based  on  experimental  data 
which  already  indirectly  include  the  coriolis  component  in  the  observations 
used.  There  is  an  exception  in  inertia  currents  where,  theoretically,  the 


I 


coriolis  force  is  the  only  acting  orce.  Coriolis  force  is  othervise  a  modi-  j 

xying  factor  and  will  be  briefly  c  nsidered  in  tho  next  chapter.  j 
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3.  EVALUATION  OF  MODIFYING  FACTORS 

The  modifying  factors  of  surface  currents  are  here  defined  as  the  factors 
vhich  do  not  directly  cause  a  current  but  affect  its  speed  and  direction  both 
in  the  horizontal  and  vertical.  In  addition,  this  chapter  considers  the  change 
of  current  speed  with  depth  in  general  and  considers  special  relationships  be¬ 
tween  surface  and  subsurface  ("deep")  currents. 

3.1  Change  of  Current  Velocity  with  Depth 

Any  current  speed  and  direction  change  with  depth  may  be  related  to  inter¬ 
nal  friction,  in  connection  with  different  currents  in  deeper  water,  bottom 
friction  or  viscosity.  The  change  of  current  with  depth  can  also  be  related 
to  such  features  as  density  gradients  with  depth  (stability),  depth  of  the 
mixed  layer  and  interactions  between  inertia  and  driving  forces  at  the  sur¬ 
face.  No  exact  model  of  current  speed  change  with  depth  exists. 

Ekman  [25]  derived  a  theory  for  the  change  of  current  speed  and  direction 
with  depth  at  the  beginning  of  the  century.  This  theory  was  later  modified  by 
Rossby  and  Montgomery,  who  related  the  deflection  to  latitude.  Available  mea¬ 
surements  show  that  Etansn's  theoretical  consideration  and  the  Rossby -Montgomery 
modification  do  not  verify  in  the  ocean,  a  fact  recognized  by  Ekman  himself. 
However,  no  available  measurement  series  permits  construction  of  a  generally 
valid  theoretical  or  empirical  model  of  the  current  structure  with  depth,  or  a 
direct  analysis  of  the  factors  involved,  because  of  the  shortness  of  the  mea¬ 
surement  series,  lack  of  secondary  data  and  great  variability  of  existing  data 
in  space  and  time. 

On  the  basis  of  theoretical  considerations,  the  various  components  of  sur¬ 
face  currents  differ  considerably  with  depth.  The  pure  tidal  current  should 
have  uniform  motion  from  the  surface  down  (except  very  close  to  the  bottom). 

The  depth  change  of  the  wind  current  without  friction  should  be  a  straight  line 
with  only  a  slight  decrease  of  speed  with  depth.  In  an  inclined  wind  drift  with 
friction,  the  depth-speed  graph  is  a  parabola.  The  depth  change  of  mass  trans¬ 
port  by  waves  is  exponential.  The  permanent  flow  can  be  considered  to  have 
uniform  speed  with  depth  in  the  mixed  layer  but  a  relatively  large  speed  gra¬ 
dient  at  the  pycnocline. 

Available  current  measurement  series  indicate  that  the  current  velocity 
change  with  depth  in  the  surface  mixed  layer  is  variable  in  space  and  time. 

The  greatest  changes  with  depth  occur  near  the  surface,  near  the  bottom  and  at 
the  pycnocline.  Some  of  these  changes  have  been  summarized  by  Stevenson  [105]  , 
His  data  have  been  utilized  with  the  following  assumptions  and  considerations 
in  the  construction  of  figure  3.1:  (a)  it  is  assumed  that  wind  currents  reach 
the  depth  of  the  thermocline;  (b)  the  current  speed  falls  to  about  30%  of  its 
surface  value  at  about  3  m.  above  the  bottom  or  above  a  sharp  pycnocline  (aver¬ 
age  value  of  current  at  3  meters  from  the  surface  being  100%).  The  upper 
curved  pazi;  of  the  depth-speed  graph  (fig.  3.1)  Is  caused  by  mass  transport  by 
waves . 
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Figure  3.1  Approximate  “Average”  Change  of  Wind  Current  Speed  with  Depth  in  Shallow  Water 
(after  Stevenson  (105]  .modified) 

3.1  Coastal  and  Bathymetric  Influences 

The  currents  near  the  coast  can' be  classified  Into  two  systems:  (a)  coast¬ 
al  current  systems,  consisting  of  relatively  uniform  drift,  approximately  paral¬ 
lel  to  the  shore,  composed  ?  tidal,  wind  and  gradient  currents,  and  (b)  near¬ 
shore  current  systems,  consisting  of  coast-ward  currents  at  the  surface  ("wave 
current" ) ,  a  longshore  current  and  seaward  flow  along  the  bottom  and  in  the 
current  rips. 

The  changing  depth  on  the  continental  slope  and  shelf  also  influences  cur¬ 
rent  speed  and  direction.  Fleming  and  Hagerty  [28]  found  that,  in  general,  the 
currents  have  a  tendency  to  flow  along  depth  contours  where  the  depth  changes 
more  suddenly  rather  than  loilowing  the  configuration  of  the  coastline.  Usual¬ 
ly  the  strong  permanent  currents  flew  along  the  continental  slope. 

Currents  flowing  along  the  coast  with  a  slight  "offshore  angle"  cause  up- 
welling  and  Intensive  mixing  off  the  slopes.  Eddying  and  upwelllng  are  also 
caused  by  steady  currents  along  the  coast  if  the  direction  of  the  coast  changes 
suddenly  (see  fig.  3.2).  A  complicated  configuration  of  the  coast  also  com¬ 
plicates  current  patterns  and  requires  detailed  local  evaluation.  Islands  in  a 
steady  current  cause  eddies  on  the  leeward  side  and  modify  the  current  pattern 
at  a  considerable  distance  from  the  coast.  In  semi closed  bays,  theire  Is  usually 
a  "conformal  current"  along  the  coast,  very  similar  to  the  currents  observed  in 
lakes.  The  tidal  currents  are  usually  the  dominant  currents  on  the  continental 
slope.  Upon  these  tidal  currents  there  is  a  wind  current  component  and  near  the 
estuaries  a  haline  gradient  component  caused  by  runoff.  True  permanent  flew 
can,  in  many  cases,  be  neglected  on  the  continental  shelf. 
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,2  Examples  of  Eddy  Formation  Around  Headlands  (Generali  zed  from  Laevastu,  Avery,  and 
Cox  [61]  ) 
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The  near-shore  current  system  Is  Illustrated  In  figure  3.3.  Longshore 
currents  can  be  approximately  predicted  along  straight  coasts  (U.S.  Navy  Weath¬ 
er  Research  Facility  [118]  j  Hunt  [45] ;  Putnam,  Munk  and  Traylor  [92] .  However, 
on  irregular  coasts,  they  depend  on  topography  and  are  greatly  influenced  by 
sand  bars,  channels,  etc.  Figure  3.^  gives  a  nomograph  for  predicting  long¬ 
shore  current.  This  nomograph  is  valid  only  for  straight,  uniformly  sloping 
beaches  with  parallel  contours. 

Another  formula  has  been  proposed  for  prediction  of  longshore  currents 
caused  by  mass  transport  by  waves: 

V  =  gmT  sin  2  0b  (25) 

vhere  V  is  longshore  current  velocity,  g  ,  acceleration  of  gravity,  m  ,  beach 
slope,  T  ,  wave  period  and  <?b ,  breaker  crest  angle.  No  proper  evaluation  of 
this  formula  has  been  made. 

No  methods  exist  for  predicting  longshore  currents  on  Irregular  and 
steep  beaches,  lhe  local  knowledge  and  experiences  from  similar  locations  and 
conditions  forms  the  basis  for  any  attempt  to  forecast  near-shore  currents  (tid¬ 
ed  component  excluded) . 
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3,3  Current  Boundaries,  Eddies  and  Deep  Currents 

The  oceanic  current  boundaries  and  large-scale  eddies  associated  with 
these  boundaries  have  been  described  by  several  investigators  -  Uda[114], 
Fuglister  and  Worthington  [29]  and  others.  These  current  features  must  be 
taken  into  consideration  in  synoptic  current  analysis  and  prediction. 

The  types  of  current  boundaries  are:  (a)  dynamic  (divergences  or  conver¬ 
gences  of  more  or  less  permanent  current  systems);  (b)  topographic  (caused  by 
the  topography  of  the  bottom  or  coast);  end  (c)  combined  eddy  systems.  The 
moat  usual  locations  of  the  current  boundaries  are:  (a)  near  meteorological 
fronts;  (b)  on  continental  slopes;  and  (c)  around  islands,  capes  and  banks 
(local  boundaries).  Seme  minor  surface  convergences  at  continental  slopes  can 
also  be  caused  by  internal  vaves.  The  bt ondaries  are  frequently  marked  oy; 
current  rips,  accumulation  of  flotsam  (surface  slicks),  modified  waves,  roar¬ 
ing  noises  and  fog.  Usually  there  Is  also  a  change  of  water  color  at  the 
boundaries.  The  surface  temperature  gradient  on  the  boundaries  can  be  0.5*  - 
2#  per  10  miles  but  may  be  even  larger  for  major  currents.  McLellan  [74]  de¬ 
scribes  temperature  changes  off  Nova  Scotia  of  3"  C,  per  meter  of  depth  and 
lateral  changes  of  1.5°  C.  per  100  meters. 

Convergences  of  tidal  current  (tidal  rips)  can  be  rather  extensive  in 
some  areas.  Such  tidal  rips  are  often  described  by  seamen;  one  such  descrip¬ 
tion  taken  from  the  Marine  Observer  (Vol.  27,  No.  176,  p.  203,  1957)  is  repro¬ 
duced  here: 

"On  the  7th  of  November  1956,  between  0900-0930  S.M.T.,  the  vessel  sailed 
through  a  series  of  tide  rips.  These  rips  extended  for  about  10  miles  on 
either  side  of  the  vessel  in  a  058®  -  238*  direction  and  gradually  got  progres¬ 
sively  rougher,  the  final  line  having  the  appearance  of  a  force  6  wind  on  the 
sea.  Each  line  was  about  100  yd.  in  width,  with  a  100  yd.  stretch  of  smooth 
water  in  between  lines.  On  passing  through  the  final  rip,  the  vessel  swung 
10*  to  starboard.  The  rips  she.  -*d  clearly  on  radar,  and  the  P.P.I.  screen  had 
the  appearance  of  a  neatly  ploughed  field.  There  was  no  change  in  sea  temper¬ 
ature.  Air  temperature  79*  F.,  sea  80°.  Light  variable  winds.  Position  of 
ship:  10  miles  W.  of  Perlas  Islands." 

There  exists  a  need  for  mapping  the  more  permanent  convergences  and  diver¬ 
gences  of  the  oceans,  based  on  ships'  logs  and  on  synoptic  computation  of  con¬ 
vergences  and  divergences  of  the  current  fields. 

The  hydrographic  and  especially  the  dynamic  conditions  at  current  bounda¬ 
ries  are  usually  rather  complicated.  Seme  graphical  examples  are  given  in  fig¬ 
ures  3.5  and  3.6. 

The  movement  of  current  boundaries  cannot  be  predicted  subjectively  at  pre¬ 
sent  with  any  great  accuracy,  and  only  certain  rules  of  thumb  can  be  established 
for  estimation  of  their  movements  (Laevastu  [60]).  However,  recent  numerical 
synoptic  anlayses  provide  several  approaches  which  indicate  their  positions  and 
gradients  with  considerable  accuracy  (Hubert,  Slusser  and  Laevastu  [44]). 
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Figure  3.5  Scheme  of  I'pwelling  and  Sinking  and  Resulting  Subsurface  Thermal  Structure  at  a 
Meandering  Boundary 


The  current  boundaries,  especially  convergences,  are  associated  with  mean- 
derlngs  and  eddies  whlcn  cause  sinking  and/or  upwelllng  of  deeper  water  (figs. 
3.5  and  3.6).  The  greater  the  speed  difference  on  both  sides  of  the  boundary 
(and  also  the  more  irregular  the  bathymetry  of  the  bottom)  the  more  extensive 
is  the  meandering.  Many  of  the  large-scale  meanderers  are  stationary,  due  to 
the  topography  of  the  continental  slope.  However  their  intensity  seems  to 
fluctuate,  depending  on  changes  in  driving  forces. 

There  are  also  great  eddy  systems  which  are  caused  by  winds  (cyclones  ar.d 
anti -cyclones) .  In  general,  the  centers  of  cyclonic  eddies  (anti-clockwise) 
are  cold,  and  they  are  associated  with  upwelllng  in  the  center.  The  centers  of 
anti-cyclonic  eddies  are  warm,  and  sinking  occurs  very  slowly  In  the  center.  If 
a  cyclonic  eddy  is  cut  off  from  a  cold  water  mass  into  a  warmer  one,  slow  sink¬ 
ing  takes  places  also  in  this  eddy  (fig.  3.7). 

Several  meteorologists  and  oceanographers,  notably  Pal  men  and  Rossby,  have 
described  the  similarities  between  atmospheric  and  oceanic  structural-dynamic 
features.  Newton  [85]  also  found  that  oceanic  and  atmospheric  eddies  are  compa¬ 
rable  in  the  sense  that  in  each  case  the  maximum  subsidence  of  the  cold  mass 
took  place  at  about  half  the  rate  of  horizontal  flow  of  mass  into  the  cold-water 
or  cold-air  tongue. 

The  currents  in  deeper  layers  of  the  sea  are  not  subjects  of  this  summary. 
However,  same  of  their  characteristics  must  be  borne  in  mind  when  analyzing  and 
forecasting  surface  currents  inasmuch  as  they  Influence  each  other.  The  deep 
currents  balance  the  "water  budget"  and  bring  to  equilibrium  the  disturbances 
of  the  surface  layers.  Therefore,  the  deep  currents  may  often  be  opposite  in 
direction  to  those  at  the  surface.  The  structure  of  deep  currents  with  depth 
can  be  complicated  (as  shown  by  recent  measurements  with  neutrally  buoyant 
floats).  Furthermore,  these  measurements  Indicate  that  the  deep  currents  are 
relatively  unsteady  in  direction  and  speed,  contrary  to  earlier  beliefs. 


SECTION  C  -  0  TEMPERATURE  CHANGES  WITH 

DEPTH  AT  C  ANO  0 


A 


j  WARMER 

B 
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SECTION  A  -  B  TEMPERATURE  CHANGES  WITH 

OEPTH  AT  A  AND  B 


Figure  3,6  Structural  Models  of  Convergences  and  Divergences 


20  ■■ 


The  well-known  coriolis  force  as  a  modifying  factor  has  no+.  been  mentioned 
earlier;  it  i&  usually  indirectly  included  in  the  semi-empirical  approaches  of 
current  analyses  and  forecasting. 


v  1  Figure  3,7  Mechanism  of  Cutting  Off  Eddies  at  Current  Boundary  and  Their  Dissipation 
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4.  OBJECTIVE  COMPUTER  METHOD  FOR  ANALYSIS 
AND  FORECASTING  OF  SURFACE  CURRENTS 

The  detail  and  accuracy  of  analyses  and  forecasts  depend  greatly  on  the 
area  and  time  scales  used.  This  is  true  both  for  objective  and  subjective  ap¬ 
proaches.  Small-scale  analyses/predictions  require  one  to  program  local  know¬ 
ledge  of  influencing  factors  as  well  as  more  detail  of  the  driving  forces.  In 
this  chapter,  we  describe  briefly  the  principles  of  hemispheric  (or  oceanwid*) 
analysis  and  prediction  of  surface  currents  on  a  12-hour  synoptic  schedule. 

The  computation  of  permanent  current  is  made  with equations(4) and (5)  (sec¬ 
tion  2.1).  In  recent  programs  the  average  temperature  between  0  and  200  meters 
is  found  by  integration,  using  FNWF  standard  levels  of  subsurface  temperature 
analyses.  As  the  analysis  refers  to  surface  currents,  and  as  countercurrents 
are  known  to  exist  below  the  permanent  pycnocline,  there  is  no  need  to  inte¬ 
grate  the  temperature  field  below  200  m.  level.  Salinity  corrections  are  added 
on  a  seasonal  basis. 

The  wind  current  component  is  computed  with  equation(7)and  added  to  the 
permanent  flow.  Integrated  true  surface  winds  are  used. 

According  to  Ekman  [25]  ,  the  direction  of  the  wind  current  at  the  surface 
is  45“  to  the  right  of  the  wind  in  the  northern  hemisphere,  and  this  angle  in¬ 
creases  with  depth.  Recent  investigations  reveal  that  the  deflection  is  more 
nearly  12-20  degrees,  being  larger  and  more  irregular  at  lower  wind  speeds 
(possibly  because  of  the  increased  importance  of  other  components)  and  smaller 
at  higher  wind  speeds.  As  the  surface  wind  is  about  the  same  angle  to  the  left 
of  the  geostrophic  wind,  it  can  be  assumed  that  the  direction  of  the  wind  cur¬ 
rent  is  the  direction  of  the  geostrophic  wind. 

It  is  assumed  that  the  current  is  relatively  uniform  and  unidirectional  in 
the  turbulent  mixed  layer  down  to  the  thermocline  (or  aoout  200  meters).  The 
mass  transport  due  to  waves,  however,  modifies  this  picture  as  this  effect  de¬ 
creases  exponentially  with  depth  (Masch  [73]).  Therefore,  if,  in  equation(7i  V 
is  in  m./sec.  and  W  in  cm. /sec.,  Kzis  taken  to  be  4,8  for  surface  currcr.ta 
(ship  routing  and  drift  computations )  and  3.2  for  the  average  current  down  to 
the  thermocline  (for  convergence/divergence  computations ) .  Obviously,  there  is 
a  time  lag  between  the  change  of  the  wind  and  response  of  the  sea.  This  lag 
seems  to  be  shorter  than  previously  believed  and  is  partially  minimized  by  the 
24-hour  averaging  of  wind  speed,  (see  further,  chapter  6). 

Since  all  computations  are  carried  out  in  the  standard  Fleet  Numerical 
Weather  Facility  grid  system,  u  and  v  current  components  are  determined  at  ap¬ 
proximately  200  nautical  mile  intervals  for  all  northern  hemisphere  ocean  areas. 
From  these  components,  direction  and  total  transport  (nm./day)  fields  are  deter¬ 
mined  and  stored  on  magnetic  tape  for  later  output  in  chart  form  or  as  special 
messages  giving  the  currents  at  specified  latitude/longitude  intersections . 

Some  current  transport  charts  are  shown  in  figures  4.1  and  4.2. 

The  contours  represent  total  current  transport  in  nautical  miles  per  day. 
One  can  clearly  Extinguish  such  well-known  features  as  the  Gulf  Stream, 


•  V* 


Figure  4.1  Surface  Current  Transport  in  Nautical  Milos  Per  Day  on  1SZ  23  May  1965 


Figure  4.2  Surface  Current  Transport  in  Nautical  Miles  Per  Day  on  00Z  25  May  1985 
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Labrador  Current,  Kuroshio  and  Cyashio.  The  low-latitude,  westerly  return  flow 
(Equatorial  current)  which  results  primarily  from  the  "wind  component"  term  is 
well-defined  in  both  the  Atlantic  and  Pacific.  A  narrow  equatorial  countercur¬ 
rent  was  obtained,  as  a  result  of  the  200  meter  temperature  structure  used  in 
the  characteristic  component . 

The  charts  ere  drawn  automatically  on  an  incremental  x-y  curve  plotter. 

Ea.cn  chart  requires  approximately  one  minute  to  complete  and  is  of  sufficient 
quality  that  it  can  be  used  immediately  for  radio-facsimile  transmission. 

In  order  to  obtain  a  single  continuous  field  display  of  both  direction 
and  speed  of  the  computed  currents,  a  stream  function  (if/)  analysis  is  made,  us¬ 
ing  methods  similar  to  those  employed  by  Bedient  and  Vederoan  [7]  ,  to  present 
atmospheric  flow  in  the  tropics.  The  vorticity  of  the  current  flow  is  deter¬ 
mined  from  the  (u,  v)  component  fields  and  the  Poisson  equation. 

,  JjL  _  JjL  (26) 

d  x  d  y 

is  solved  for  if/  using  relaxation  techniques. 

The  stream  function  fields  which  correspond  to  the  current  transport  charts 
on  figures  4.1  and  4.2  are  shown  in  figures  4.3  and  4.4.  If  current  vectors  are 
plotted  on  these  charts,  there  Is  a  very  dose  correspondence,  but  slight  cross 
contour  flow  can  be  observed.  This  is  natural  because  the  derived  stream  func¬ 
tion  is  nondivergent,  while  there  is  divergence  in  the  initial  velocity  field. 

In  general,  the  divergence  is  small  in  most  places,  and  the  stream  field  pro¬ 
vides  a  good  representation  of  the  current  pattern. 

It  can  be  seen  that  the  hemispheric  grid  size  is  too  large  to  show  details 
in  areas  of  strong  currents.  Detailed  models  sire  under  development  to  remedy 
this  shortcoming  in  areas  of  relatively  dense  synoptic  weather  reports.  Some 
of  these  models  will  also  Include  inertia  effects  and  tidal  currents.  The  re¬ 
sults  of  these  numerical  analyses/foreoasts  are  used  as  inputs  in  a  number  of 
other  oceanographic  analyses /forecasts. 


*- idea 5*n*'*.' -w^t*** v  ■■ 
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5.  VERIFICATION  OF  SURFACE  CURRENT  ANALYSES 

A  synoptic  current  chart  would  be  of  little  value  if  it  could  not  be  veri¬ 
fied  and  the  computational  a c heme  tuned  as  required.  Direct  current  measure¬ 
ments  in  the  open  ocean  are  few,  and  drift  calculations  made  from  navigational 
fixes  are  frequently  inaccurate  In  weak  current  areas,  so  that  It  is  difficult 
make  a  direct  evaluation.  It  has  been  necessary,  therefore,  to  resort  to  in¬ 
direct  means  which  are  susceptible  to  verification  on  a  synoptic  basis. 

Sea-durface  temperature  (SST)  la  the  only  oceanographic  element  which  per¬ 
mits  a  reasonably  complete  synoptic  analysis  on  a  hemispheric  scale.  Such  anal¬ 
yses  are  made  twice  dally  at  Fleet  Numerical  Weather  Facility,  Monterey  (Wolff 
(1281 ) ,  and  their  resolution  la  such  that  SST  changes  car.  be  determined  for 
periods  of  24,  48  hours,  etc.  From  these  changes  are  subtracted  the  local 
changes  computed  from  air/sea  heat  exchange  equations.  If  the  remainder  cerv> 
latea  well  with  the  advective  change  indicated  by  W  -V  SST  the  computed 
currents  can  be  assuaed  to  be  reasonably  correct.  *rz 

Same  experimental  verification  has  been  carried  out  at  Fleet  Numerical 
Weather  Facility.  The  results  have  been  very  encouraging.  These  first  compu¬ 
tations  indicate  that  the  advectlonal  change  of  SST  is,  over  major  parts  of  the 
ocean,  <0.1Cy{24h)  .  In  only  a  few  areas  does  the  change  rise  above  0.3°C.X24h) 
higher  values  are  confined  to  sharp  boundaries  of  major  currents.  Same  dis¬ 
crepancies  between  analyzed  and  computed  quantities  have  also  been  ascertained; 
these  have  led  to  plans  of  shortening  the  SST  analysis  period  to  decrease  the 
apparent  lag.  Both  SST  adveetiom  and  heat  exchange  will  be  used  In  SST  fore¬ 
casts.  Furthermore,  100  nautical  miles  and  smaller  mesh  lengths  will  be  used 
in  the  future  for  current  analyses  to  make  them  directly  compatible  to  the  SST 
analyses  which  are  made  on  a  100  nautical  mile  grid. 


6.  THE  RESPONSE  TIME  OF  SURFACE  LAYERS  TO  TRANSFER 
OF  MECHANICAL  ENERGY 


As  the  surface  layer  (which  is  set  in  motion  by  wind)  has  considerable 
Inertia,  it  can  be  expected  that  the  response  of  the  sea  to  the  change  of  driv¬ 
ing  force  is  gradual,  and  that  there  Id  a  tine  lag  between  the  change  of  wind 
speed  and  direction  *»n<i  the  response  of  the  surface  current  to  tnis  change. 

A  lag  also  exists  in  the  tine  for  growth  of  a  fully  arisen  sea  in  relation  to 
fetch  length  and  duration  of  the  wind.  S3ie  available  data  do  not  allow  an  ex¬ 
act  numerical  evaluation  of  the  response  time  but  only  &  qualitative  analysis. 
These  data  indicate  that  the  response  time  for  currents  varies  with  the  wind 
speed  (especially  with  the  ratio  between  previous  and  new  (actual]  wind  speed), 
the  depth  of  the  water,  distance  from  and  configuration  of  the  coast.  Krause 
(personal  communication)  has  also  found  a  rapid  response  in  the  relatively  shal¬ 
low  Kiel  Bight. 


A  few  additional  observations  fretn  areas  with  slow  permanent  flow  indicate 
that  the  complete  response  time  with  moderate  winds  is  about  tnree  hours.  With 
strong  winds  end  in  areas  of  strong  currents,  the  response  is  gradual,  and  the 
establishment  of  a  steady  state  takes  a  longer  time  (up  to  la  hours). 


It  should  be  underlined  here  that  the  current  direction  and  speed  In  a 
given  position  in  the  open  ocean  might  not  be  directly  determined  by  winds  pre¬ 
vailing  at  this  position  but  by  stronger,  perhaps  more  persistent  and  larger 
wind  fields  some  distance  away  from  the  given  position. 


Hafner  [33 j  offers  the  io±iowing  empirical  formula  for  estimation  of  com¬ 
plete  response  time  in  shallow  water. 


Response  time  (h) 


Fetch  length  (km) 

Vg  x  depth  of  water  (m) 


This  formula  would  merit  testing  in  deep  water.  If  we  assume  the  depth  to  be  the 
depth  of  the  surface  mixed  layer. 


It  can  be  expected  that  verification  attempts  with  numerical  analyses/fore¬ 
casts  will  shed  more  quantitative  light  on  the  response  time. 
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7.  USE  OF  CURRENT  FORECASTS  IN  NAVAL  PROBLEMS; 

NAVIGATION  AND  FISHERIES 

The  analysis  and  forecasting  of  surface  currents  on  a  synoptic  schedule 
were  initiated  at  Fleet  Numerical  Weather  Facility  primarily  because  of  their 
importance  in  the  following  problems; 

(a)  Quantitative  computation  of  divergence  and  convergence  and  the  accom¬ 
panying  up  and  down  movement  of  thermocline  depth; 

(b)  Computation  of  the  decay  of  transients  where  one  of  the  prime  factors 
is  current; 

(c)  Forecasting  the  advective  part  of  sea  surface  temperature  changes; 

(i)  Construction  of  detailed  dynamic  models  of  thermal  structure  in  areas 

of  strong  currents. 

Current  analyses  and  forecasts  also  have  other  applications.  At  present, 
navigators  use  Information  on  monthly  average  currents  available  in  atlases 
and  information  on  tidal  currents  given  on  sea  charts, in  tidal  current  tables 
and  in  pilot  books.  In  many  instances,  however,  it  seems  to  be  economically 
profitable  to  use  prognoses  of  actual  currents,  especially  in  sea  and  weather 
routing  of  ships.  This  is  now  being  done  empirically  and  indirectly. 

At  present,  allowance  is  seldom  made  in  the  routing  practice  for  drift 
caused  by  currents,  end,  therefore,  some  inaccuracies  might  be  experienced.  In 
ship  routing  practice,  the  resultant  of  wind  currents  and  permanent  flow  could 
be  used.  The  influence  of  "wave  currents"  could  also  be  considered  as  one  or 
two  separate  factors:  (a)  involuntary  speed  reduction  and  (b)  resistance  by 
wave  currents.  Direct  use  of  information  on  currents  is  made  in  rescue  opera¬ 
tions,  in  estimations  of  drifts  of  polluted  or  otherwise  contaminated  waters 
and  in  the  estimation  of  ice  flow. 

In  almost  all  types  of  fisheries  work  there  is  a  need  to  know  the  actual 
currents  and  to  adjust  fishing  operations  accordingly.  The  need  for  and  the 
use  of  currents  in  fisheries  has  been  summarized  by  Hela  and  L&evastu  [40]  . 
Current  boundaries  are  extremely  important  to  fisheries.  At  divergences  the 
deeper,  nutrient-rich  water  is  brought  into  the  surface  layers  where  it  causes 
a  higher  production  of  organic  matter  and  an  aggregation  of  fish;  convergences 
cause  (dynamically)  a  concentration  of  zooplankton  and  are  accompanied  also  by 
a  concentration  of  fish.  Pic  tori  ally,  one  can  say  that  there  is  an  accumula¬ 
tion  of  everything  from  plankton  to  fishermen  on  a  convergence. 

Numerical  programs  for  analysis  and  prediction  of  current  and  water  type 
boundaries  at  the  surface  are  under  development  and  testing  at  FNWF.  These 

utilize  SsT  gradients  zc  well  p-  surface  ccfttut  ,*i-*sdients . 

Past  knowledge  of  actual  surface  currents  has  been  shaky  indeed;  this  fact 
has  prevented  their  effective  economical  use  in  navigation  and  fisheries  work. 
The  development  of  good,  synoptic  analyses  and  forecasts  of  ocean  currents  would 
benefit  many. 


FRECSDItC  PASS  BLANK,  HOT  FIXMEDj 
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APPENDIX 


SYMBOLS  AND  UNITS  USED 
A  Amplitude  of  the  tide 

Ac  Current  component  caused  by  the  change  of  atmospheric  pressure 
B  Breadth  of  current 

C 

C  Wave  speed 

AD  Dynamic  height  difference  (in  dynamic  cm.) 
d  Depth  of  water 
f  Corioxia  parameter 

G  Permanent  flow  (gradient  current  thermohaline  and/or  characteristic 
c  current  large-scale,  wind-driven  ) 

g  Acceleration  of  gravity 

AHg  oeostrophic  potential  (/3  -  «a  ph) 

AH ^  Longitudinal  height  difference  (of  sea  level) 

AHt  Transversal  height  difference  (of  sea  level) 

AHm  Slope  of  sea  surface  (cm./lOO  km.) 

H1/sWave  height  (usually  significant  height) 

h  L  Thickness  of  top  layer 

hs  Thickness  of  lower  (bottom  layer) 

I  Velocity  and  directional  component  of  surface  current  caused  by  influ¬ 
encing  factors 

i  Inclination  of  the  sea  surface 

Kj  K2  Constant,  factor 

jK  Unit  vector 

L,/3  Wave  length 

Distance  between  stations; 

1  Distance  in  nautical  miles; 

Length  of  basin 

A-l 


m  Beach  slope 
p  Pressure 

F,  Periodic  port!  on  of  inertia  current 
P,  Periodic  portion  of  tidal  current 
r  Radius  of  inertia 
T  Wave  period 

T  Mean  \  snperature  above  zero  current 
T  Period  of  inertia 

p 

T  Sea-surface  temperature 
t,t  Period  of  internal  waves 
V  Wind  speed  (m./sec.) 

Vt  Velocity  of  internal  waves 
W  Current  speed 
W  Wind  currerkt 

c 

Longitudinal  velocity  of  current 
Wm  Mass  transport  velocity  ("wave  current") 

VYpls,  Velocity  of  surface  current,  with  given  direction  (place,  time  and  depth 
of  the  actual  current  being  specified) 

Wj-,  Wa  Relative  current  speeds  in  two  differs  art  layers 

x  Friction  constant  (bottom  friction;  O.OQ25  -  0.3  usually  taken  0.3) 

As  Depth  of  aero  current 

Specific  volume  anomaly  (constant  along  °t  surface) 
sin  «  Slope  of  the  bottom 

/3  Anomaly  of  dynamic  height 
/3  Deflection  of  the  wind  current  (in  p) 
tan  y  Downward  inclination  of  the  pycnoclme 


A-2 


l 

( 


<P  Average  geographical  latitude 
8b  Breaker  crest  an gfte 

a)  Angular  velocity  of  the  earth  (7*29  10  rad./aec.) 

Densities  of  two  layers  (of  top  and  lower  layers) 
A (>  Density  difference  atrtnTerface 
r  Shearing  stress  of  wind 
<7(  Sea  water  density  temperature 
\]j  Stream  function 
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Results  of  large  scale,  synoptic  current  calculations  I  Results  ofiarge  scale.  :  .ynoptlc  current  calculation 

are  shown.  j  are  shown. 
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